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ABSTRACT 

We discuss a method to constrain the intrinsic shapes of galaxy clusters by combining X-Ray and 
Sunyaev-Zeldovich observations. The method is applied to a sample of 25 X-Ray selected clusters, 
with measured Sunyaev-Zeldovich temperature decrements. The sample turns out to be slightly biased, 
with strongly elongated clusters preferentially aligned along the line of sight. This result demonstrates 
that X-Ray selected cluster samples may be affected by morphological and orientation effects even if 
a relatively high threshold signal-to-noise ratio is used to select the sample. A large majority of the 
clusters in our sample exhibit a marked triaxial structure; the spherical hypothesis is strongly rejected 
for most sample members. Cooling flow clusters do not show preferentially regular morphologies. We 
also show that identification of multiple gravitationally-lensed images, together with measurements of 
the Sunyaev-Zeldovich effect and X-Ray surface brightness, can provide a simultaneous determination 
of the three-dimensional structure of a cluster, of the Hubble constant, and the cosmological energy 
density parameters. 

Subject headings: Galaxies: clusters: general - X-Rays: galaxies: clusters - cosmology: observations 
- distance scale - gravitational lensing - cosmic microwave background 



1. INTRODUCTION 

The intrinsic, three-dimensional (hereafter 3-D) shape 
of clusters of galaxies is an important cosmological probe. 
The structure of galaxy clusters is sensitive to the mass 
density in the universe, so knowledge of this structure 
can help in discriminating between different cosmological 
models. It has long been clear that the formation epoch 
of galaxy clusters strongly depe nds on the matter den - 
sity parameter of the universe l)Richstone et all 11992ft . 
The growth of structure in a high-mattcr-density uni- 
verse is expected to continue to the present day, whereas 
in a low density universe the fraction of recently formed 
clusters, which are more likely to have substructure, is 
lower. Therefore, a sub-critical value of the density pa- 
rameter f^Mo favors clusters with steeper density profiles 
and rounder isodensity contours. Less dramatically, a 
cosmological constant also delays the formation epoch 
of c lusters, favoring t he presence of structural irregular- 
ity l)Suwa et al.ll2003ft . 

An accurate knowledge of intrinsic cluster shape is also 
required to constrain structure formation models via ob- 
servations of clusters. The asphericity of dark halos af- 
fects the inferred central mass density of clusters, the 
predicted frequency of gravitational arcs, nonlinear clus- 
tering (especially high-order clus tering statistics) a nd dy- 
namics of galactic satellites (see lJing fc Sutol i|2002ft and 
references therein). 

Asphericity in the gas density distribution of clusters of 
galaxies is crucial in modeling X-Ray morphologies and 
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in usi ng clust ers as cosmolog ical tools. l)Inagaki et alJ 
119951 ICooravl 119981 ISulkanenl E299) . Assumed cluster 
shape strongly affects absolute distances obtained from 
X-Ray/Sunyaev-Zeldovich (SZ) measurements, as well 
as relativ e distances obtained from baryon fraction con- 
straints l|Allen et alJ 120041: ICooravl 11998ft . Finally, all 
cluster mass measurements derived from X-Ray and dy- 
namical observations are sensitive to the assumptions 
about cluster symmetry. 

Of course, only the two-dimensional (2-D) projected 
properties of clusters can be observed. The question 
of how to deproject observed images is a well-posed 
inversi on problem that has been studied by m any au- 
thors (Ejc3[1971 E^de^[l99l Eeblinik3|2003). Since 
information is lost in the process of projection it is in 
general impossible to derive the intrinsic 3-D shape of 
an astronomical object from a single observation. To 
some extent, however, one can overcome this degener- 
acy by combin ing observations in differen t wavelengths. 
For example, iZaroubi et alJ l)1998l 12001ft introduced a 
model-independent method of image deprojection. This 
inversion method uses X-Ray, radio and weak lensing 
maps to infer the underlyin g 3-D structure f or an axi- 
ally symmetric distribution. iReblinskvl l)2000ft proposed 
a parameter-free algorithm for the deprojection of ob- 
served two dimensional cluster images, again using weak 
lensing, X-Ray surface brightness and SZ imaging. The 
3-D gravitational potential was assumed to be axially 
symmetric and the inclination angle was required as an 
input parameter. Strategies for determining the o rien- 
tation have been also discussed. iDore et alJ 1)200 fj) pro- 
posed a method that, with a perturbative approach and 
with the aid of SZ and weak lensing data, could predict 
the cluster X-Ray emissivity without resolving the full 
3-D structure of the cluster. The degeneracy between 
the distance to galaxy clusters and the elongation of the 
cluster alon g the line of sight (l.o.s.) was thoroughly 
discussed bv lFox &: Pen! <)2002ft . They introduced a spe- 
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cific method for finding the intrinsic 3-D shape of triaxial 
cluster and, at the same time, measuring the distance to 
the cluster corrected for asphericity, so prov iding an un- 
biased estimate of the Hubble constant Hq. iLee fc Sutol 
( 2004) recently proposed a theoretical method to recon- 
struct the shape of triaxial dark matter halos using X- 
Ray and SZ data. The Hubble constant and the projec- 
tion angle of one principal axis of the cluster on the plane 
of the sky being independently known, they constructed 
a numerical algorithm to determine t he halo eccentrici- 
ties and orientation. H owever, neither IFox fc Pen! (2002) 
nor ILee fc Sutol (|2004fl apply their method to real data. 
In this paper we focus on X-Ray surface brightness ob- 
servations and SZ temperature decrement measurements. 
We show how the intrinsic 3-D shape of a cluster of galax- 
ies can be determined through joint analyses of these 
data, given an assumed cosmology. We constrain the tri- 
axial structure of a sample of observed clusters of galax- 
ies with measured X-Ray and SZ maps. To break the 
degeneracy between shape and cosmology, we adopt cos- 
mological parameters which have been relatively well- 
determined from measurements of the cosmic microwave 
background (CMB) anisotropy, Type la supernovae and 
the spatial distribution of galaxies. We also show how, 
if multiply- imaging gravitational lens systems are ob- 
served, a joint analysis of strong lensing, X-Rays and 
SZ data allows a determination of both the 3-D shape of 
a cluster and the geometrical properties of the universe. 
The paper is organized as follows. The basic dependen- 
cies of cluster X-Ray emission and the SZE on geometry 
are reviewed in § In § [21 we show how to reconstruct 
the 3-D cluster structure from these data, presuming cos- 
mological parameters to be known. In passing we note 
how the addition of suitable strong gravitational lensing 
data can constrain the cosmological parameters as well, 
although we do not impose lensing constraints in this pa- 
per. We then turn to face the data. Our cluster sample is 
introduced in § and in § we present 2-D X-Ray sur- 
face brightness parameters for each sample member. The 
triaxial structure of the clusters is then estimated and an- 
alyzed in § EI § d is devoted to a summary and discussion 
of the results. In Appendix 1X1 we provide details on the 
triaxial ellipsoidal /3-model, used to describe the intra- 
cluster gas distribution, while Appendix IBl is devoted to 
a discussion of the consequences of our assumption of 
clusters being triaxial ellipsoids aligned along the line 
of sight. In Appendix ICl the identifications of multiple 
sets of images of background galaxies in strong lensing 
events is discussed. Throughout this paper, unless oth- 
erwise stated, we quote errors at the 68.3% confidence 
level. 

2. MULTI- WAVELENGTH APPROACH 

In this section, we summarize the relationships be- 
tween SZ and X-Ray observables, on the one hand, and 
cluster shape and distance on the other. 

2.1. The Sunyaev-Zeldovich Effect 

The gravitational potential wells of galaxy clusters 
contain plasma at temperatures of about k^T e w 8- 
10 keV. CMB photons that pass through a cluster 
interact with the energetic electrons of its hot intra- 
cluster medium (ICM) through inverse Compton scatter- 
ing, with a probability r ~ 0.01. This interaction causes 



a small distortion in the CMB sp ectrum, known as the 
Sunyaev-Zeldovich eff ect (SZE) l)Sunvaev fc ZeldovichI 
119701 birkinshaw 1999), which is proportional to the elec- 
tron pressure integrated along the l.o.s., i.e. to the first 
power of the plasma density. The measured temperature 
decrement ATsz of the CMB is given by: 



AT SZ , .<nk B 



(1) 



Tomb m c u- j l o s . 

where T c is the temperature of the ICM, fee the Boltz- 
mann constant, Tcmb = 2.728°K is the temperature of 
the CMB, or the Thompson cross section, m e the elec- 
tron mass, c the speed of light in vacuum and f(u,T e ) 
accounts for frequency shift and relativistic corrections. 
If we assume that the ICM is described by an i sothe rmal 
triaxial /3-model distribution, substituting Eq. (|A10|) into 
with m — 1, we obtain: 

/ „9 , 9 V2-3/3/2 
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where ATo is the central temperature decrement which 
includes all the physical constants and the terms result- 
ing from the l.o.s. integration 
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with: 
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D c is the angular diameter distance to the cluster, 
9i = Xi t0 hs/D c is the projected angular position (on the 
plane of the sky) of the intrinsic orthogonal coordinate 
£i,obs, h is a fun ction of the cluster shape and orien- 
tation (Eq. IA11|) . e pro j is the axial ratio of the major 
to the minor axes of the observed projected isophotcs 
and c ,proj the projection on the plane of the sky ( p.o.s.) 
of the intrinsic angular core radius (Eq. IA12|) . In 
a Friedmann-Lemaitre-Robertson- Walker universe filled 
with pressure-less matter and with a cosmological con- 
stant, the angular diameter distance between an observer 
at a redshift z& and a source at z c is: 

i i „. f r c \n K0 
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(5) 



= Vn M o(i + z) 3 + n A0 + n K0 (i + z) 2 

where Hq, £Imo and f^Ao are the Hubble parameter, the 
normalized energy density of pressure-less matter and 
the reduced cosmological constant at z = 0, respectively. 
Qko is given by f^Ko = 1 — f^Afo - ^ao, and Sinn is defined 
as being sinh when Qko > 0, sin when f2xo < 0; and as 
the identity when f^xo — 0. A more general expression of 
the angular diameter distance, also accounting for dark 
energy a nd inhomogeneity in matt er distribution, can be 
found in lSereno et all pOOll 12002]) . 
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2.2. X-Ray Surface Brightness 

Cluster X-Ray emission is due to bremsstrahlung and 
line radiation resulting from electron-ion collisions; the 
X-Ray surface brightness Sx is proportional to the inte- 
gral along the l.o.s. of the square of the electron density: 

Sx = L, nlAeHdl (6) 

where A e # is the X-Ray cooling function of th e ICM in 
the cluster rest frame. Substituting Eq. l{A"lT)j) into © 
with m = 2, we get: 



?x = Sxo 



,2, p 2 fl2\ 1/2-3/3 
1 1 "T e V ro f2 \ 



9 2 . 

c.proj 



(7) 



where the central surface brightness Sxo reads: 



SxO = A r-„ . ~ n e0 TTTa \ l 9{P) (8) 



4^(1 + ^ c ) 4 e ° ^3/4 v , j>rii 

fi is the molecular weight given by: /i^ = p/nim p . 

3. COMBINING HETEROGENEOUS DATA SETS 

Here we discuss how 2-D SZE and X-Ray maps of a 
cluster can be used to constrain its 3-D shape. We fol- 
low a parametric approach. We model the cluster us- 
ing an isothermal, triaxial [3 profile, and adopt a concor- 
dance model for the cosmological distance-redshift rela- 
tionships. Details of the cluster model are given in Ap- 
pendix El 

This model has long been used to describe the elec- 
tron distributions of galaxy clusters. It was originally 
introduced sp ecifically for dynamically re l axed, isother- 
mal clusters ijCavaliere fc Fusco-Femianol I1978D . but it 
was then observed to fit the X-Ray emission of most 
galaxy clusters reasonably well. A serious drawback of 
this model is that electron density profiles with extreme 
axial ratios lead either to unlikely total mass density dis- 
tributions, i.e. dumbbell shaped clusters, or to regions 
with unphysical (negative) density. Nevertheless, we be- 
lieve its extreme versatility makes it a useful tool for our 
purposes. 

For an ellipsoidal distribution, the 3-D shape of a clus- 
ter is described by two axis ratios, e\ and e2, and the 
orientation of the cluster is fixed by th ree Eule r angles, 
#Eu,0Eu and Veu- As shown in Eqs. (|A1I IC2|I . in our 
model the density profile of a cluster is characterized by 
three additional parameters: the central density n e o, the 
slope (3 and a core radius, r c ^. Under the hypothesis 
of isothermal ICM, a single value, T c , characterizes the 
temperature profile of the cluster. In all, nine parame- 
ters describe the cluster. 

As discussed in the previous section, the cosmological 
dependence of the model enters through the luminosity- 
rcdshift relationship. For a flat model universe, this re- 
lationship is in turn determined by two parameters: the 
Hubble constant, Hq, and the matter density Qmo- 
A projected axis ratio, e pro j, and an orientation angle, i/j, 
characterize a family of ellipses in the p.o.s. derived from 
the 2-D projection of 3-D ellipsoids. By fitting an ellip- 
tical profile to the X-Ray and/or SZE data, these two 
parameters can (in principle) be measured. Two other 
observables, the slope /3 of the profile and the projected 
core radius # c ,proj can also be determined from data. Two 



independent ge ometrical constraints relate 2-D and 3-D 
quantities (Eqs. EH EH - 

So far we have discussed only quantities derivable from 
spatial distributions. Besides these, the cluster central 
electron density and the temperature of the ICM can 
be inferred from X-Ray observations with sufficient en- 
ergy resolution. The observed values of the central tem- 
perature decrement, ATq in Eq. and of the central 
surface brightness, Sxo m Eq. JSJ, provide two further 
constraints. If some assumption is made on the orienta- 
tion of the cluster, with eight independent equations and 
eight unknown physical parameters a full determination 
of the cluster shape can be obtained. If a rotational ellip- 
soidal morphology is chosen, a lower number of param- 
eters is needed to describe the three-dimensional shape 
of clusters; in this case no additional assumption on the 
inclination is required allowing a full determination of 
the cluster shape and orientati on. This case is tre ated in 

\l .1 12004ft . 



details in a subsequent paper 

3.1. Adding Strong Lensing Data 

We wish to point that, even though we do not do so 
in this work, strong lensing data can be combined with 
the X-Ray and SZE observations to break the degeneracy 
between the intrinsic shape of the lensing cluster and the 
cosmological parameters. If this were done, one could 
obtain simultaneous constraints on the cluster parame- 
ters and on the cosmology. 

In particular, each set of strong gravitational images 
identified in a cluster strongly constrains the mass dis- 
tribution of the lens. The convergence k depends on the 
cosmology only through the ratio of distances D cs /D s . 
Therefore k depends only the cosmological density pa- 
rameters Qi, and not on the Hubble constant Hq. The 
value of ko changes according to the redshift z s of the 
lensed source. Image systems produced by sources at 
different redshifts probe independent values of the ratio 
D cs /D s . 

Each image system provides a cons traint on central value 
of the convergence fco^s) (Ea. IC7|I . In turn, if both the 
positions and the redshift of a multiple image system 
arc know n, each measured value of ko provides, through 
Eq. I|C7|I . a further independent constraint on the cos- 
mological energy densities. Each multiply- imaged source 
provides an independent constraint which relates the cos- 
mological parameters fli to the 3-D shape and orientation 
of the cluster (/i 3 / 4 / v / eTe2). With a sufficient number 
of image systems, then, a measure of both the intrinsic 
shape and orientation of the cluster and a simultane- 
ous estimate of all cosmological parameters involved can 
therefore be performed. 

3.2. Angular Diameter Distances from X-Ray and SZE 
Observations for Triaxial Clusters 

It is of course well known that the angular diameter 
distance to a spherically symmetric cluster can be be in- 
ferred from microwave decrement and X-Ray data. The 
angular diameter distance enters the SZE and the X-Ray 
emission through a characteristic length-scale of the clus- 
ter along the l.o.s. SZE and X-Ray emission depend dif- 
ferently on the density of ICM, and therefore also on the 
assumed cosmology. A joint analysis of SZE measure- 
ments and X-Ray imaging observations, together with 
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the assumption of spheri cal symmetry, thus can yield 
the distance to the cluster (Birk inshawll999t lR eese et all 
2002). Specifically, one can solve Eqs. 10 and © for the 
angular diameter distance D c , by eliminating n C Q. 
More generally, for a triaxial cluster the inferred angular 
diameter distance takes the form: 



D = D \ ,proj h 1/2 

— ^c| Exp „ H 
"c3 



cIexd 1 ei e 2 



1/2 



(9) 



where -D c |e X p ^ s an experimental quantity given by: 
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Sxo \kBT e aJ g(/3/2) 2 6» CiProj 

A-eHO t l e/ r l H 



in*/* f(v,T e )*Tl „2 



. ri+TF (10) 

CMB °T V 1 + z c) 

Under the assumption of spherical symmetry, the 3-D 
morphology of the cluster is completely known: h — e\ — 
e2 = 1, ^Eu = <^Eu = ^Eu = and the observed major 
core radius 9 CyPTO j reduces to 9 C . Hence Eq. © becomes: 



D r = D r 



Exp 



and the cluster angular diameter distance can therefore 
be obtained directly from Eq. ljTU|) . The standard ap- 
proach in the past decades has been to take advantage 
of this possibility to estimate D c under the assumption 
of spherical symmetry, in order to constrain the un- 
derlying cosmology. Since in fact it is also true that 
D c = -Dclcosm' through Eq. |@J an estimate of H can 
be obtained if f£ m q and £l\n are known from indepen- 
dent observations (Birkinshaw 1999; iReese et al] 120021 
IMason et aTll2001|) . 

The same approach clearly cannot be applied when the 
assumption of spherical symmetry is relaxed and clus- 
ters are considered as more general triaxial systems. In 
this case an estimate of the axis ratios, shape and ori- 
entation parameters is required before D c can be com- 
puted. Conversely, if D c is known from the redshift and 
prior knowledge of the cosmology, then the X-Ray and 
SZE data can be used to constrain the 3-D morphol- 
ogy of the cluster. In this paper we will follow this lat- 
ter approach. We assume the values of fii and of Hq 
to be known; -Dclcosm can be- then determined through 
Eq. (@J. We will then use Eq. © to infer the 3-D mor- 
phology of a sample of galaxy clusters. We believe the 
cosmological distance scale is now known with sufficient 
accuracy to warrant our approach. An impressive body 
of evidence from CMB anisotropy, Type la supernovae, 
galax y clustering, large- scale structure, and the Lya for- 
est (jWang et alJ l200T]h are consistent with a the pic- 
ture of a universe with sub-critical cold dark matter 
energy density and with two-thirds of t he critical den- 
sity b eing in the form of dark energy. iTeemark et alJ 
( 2004) combine the three dimensional power spectrum 
from over 200, 000 galaxies in the Sloan Digital Sky Sur- 
vey with the first-year Wilkinson Microwav e Anisotropy 
Probe (WMAP) data (|Spergel et all 120031) to measure 
cosmological parameters. Their results are consistent 
with a flat (17k = 0) cosmological model with Hq = 
70ta km s" 1 Mpc -1 , fi M = 0.30±0.04 and with a non- 
zero cosmological constant. Thanks to the high precision 



to which cosmological parameters are known, we are able 



to constrain measurements of D, 
jects within a 5% error. 



c I Cosm 



for the sample ob- 



3.3. Cluster Elongation Along the Line of Sight 

For the remainder of this paper, we will assume that 
every cluster is triaxial, with one principal axis aligned 
along the l.o.s. (see §EJ. In Appendix iBl we show that 
the magnitude of the systematic error in inferred elonga- 
tion parameters cause by such assumption is small com- 
pared to the uncertainties arising from the observational 
errors. Such a straightforward assumption also leads to 
an extremely simple formalism to describe the resulting 
three-dimensional shape of clusters, reducing the errors 
caused by the uncertainties in the observational data. 

The assumption that the cluster is aligned along the 
l.o.s. implies: 6» Eu = <^Eu = ^Eu = 0, h=l, j — e\, k = 
0, and I = e|; (see Appendix |A"|) . We label axes so that 
major axis is parallel to the x±; then the projected axial 
ratio and core radius are: e Dro 



v 2 /vi and CiPrc 



c/vi. The angular diameter distance becomes: 



Dr = D, 



c I Cosm 



I Exp 



(11) 



We now introduce the elongation ei.o.s., defined as the 
ratio of the radius of the cluster along the l.o.s. to its 
major axis in the p.o.s., 



Cl.o.s. — 



Dr 



Exp 



D 



c I Cosm 



(z,h ,q m )' 



(12) 
(13) 



Spherical clusters have the same radius along the l.o.s. 
and in the p.o.s. and for them ei. . s . = 1. Clusters which 
are instead more or less elongated along the l.o.s. than 
in the p.o.s. will have values of ei. . s . > 1 or ei. Q . s . < 1, 
respectively. 

4. DATA SAMPLE 

We now apply the formalism described in § 13.31 to a 
sample of galaxy clusters to infer new information about 
the extent of the clusters along the l.o.s. 
We use two samples of clusters for which combined X- 
Ray and SZ analys is has already bee n reported. The 
sample discussed bv IReese et al.l i|2002|) consists of 18 X- 
Ray selected clusters with z > 0.14 and S > —15° and 
L x (0.1 - 2.4 keV) > 5xl0 44 ^ 2 erg s" 1 and for which 
high S/N detections of SZE, high-S/N X-Ray imaging 
and electron tem peratures wer e avail able. To these we 
add the sample of lMason et al.l (|2001fl , w hich contains 7 
cluste rs from X- Ray-flux-limited sample of lEbeling et alJ 
(1996). Details on the completeness of t he latter subsam- 
ple are given bv IMason fc Mversl 1)20001) . 
Basic data for our 25 clusters, including previously 
published redshift, plas ma temperature and microwave 
decre ment information ijReese et al.l 120021 IMason et alJ 
120011) are presented in Tabled 

5. X-RAY MORPHOLOGY IN TWO DIMENSIONS 

Chandra and XMM observations of clusters in the past 
few years have shown that in general clusters exhibit el- 
liptical surface brightness maps, and so cannot be spher- 
ically symmetric. In order to obtain a uniform set of 
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Table 1. Clusters in the Sample 



Cluster 


z 


fe B T„ 


A To 






(keV) 


O*0 



MS 1137.5+6625 
MS 0451.6-0305 
CI 0016+1609 
RXJ1347. 5-1145 
A 370 

MS 1358.4+6245 

A 1995 

A 611 

A 697 

A 1835 

A 2261 

A 773 

A 2163 

A 520 

A 1689 

A 665 

A 2218 

A 1413 

A 2142 

A 478 

A 1651 

A 401 

A 399 

A 2256 

A 1656 



0.784 
0.550 
0.546 
0.451 
0.374 
0.327 
0.322 
0.288 
0.282 
0.252 
0.224 
0.216 
0.202 
0.202 
0.183 
0.182 
0.171 
0.142 
0.091 
0.088 
0.084 
0.074 
0.072 
0.058 
0.023 



0.8 



5 7^ ° 
°- < -0.7 

10 4+ 10 

7 cr + 0.72 

9 3+°' 7 
« K + 0-7 

6 - 6 -o.5 

' •^°-0.42 

8 sg+o- 86 

6 - 6 -0.6 
+ 0.7 
-0.7 
s 91 +0.19 
S - 21 -0.17 

o on"l"D ■ 3 ? 

8.82_ Q 32 

y-^y_0.36 
, + 1.1 
-0.7 
, + 0.46 
-0.40 
. + 0.22 
" -0.20 
, + 0.35 
-0.31 
-+0.22 
-•-0.21 
.+0.17 
-0.16 



12.27 
8.331 
9.66j 
9.03j 
7.05j 
7.54j 

7.0 ± 0.2 

8.0 ± 0.2 
8.4 ± 0.7 
6.4 ± 0.2 

9.1 ± 0.4 
9.7 ± 0.8 
6.6 ± 0.2 



-818 
-1431 
-1242 
-3950 
-785 



+ 98 
113 
98 
105 
+ 105 
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+ 350 
350 
+ 118 
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-1410+ 160 

^^"-180 
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150 
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+ 110 
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+ 30 
-30 
+ 20 



-728 

-731 

-856 
-437 
-375 
-247 
-338 _ 
-1M±* 
-243if 9 
-30211? 




Fig. 1. — Comparison of central surface brightness (top), slope 
j3 (middle) and core radius (bottom) obtained fitting circular and 
elliptical /3-models to the X-Ray surface brightness maps of sample 
clusters. In the bottom panel the gray symbols show core radii 
corrected for the ellipticity as in Eq. 1151 



Note. — Clusters in the sample; their rcdshift. gas temperature and 
central temperature decrement. 



X-Ray observables for our sample objects, we have re- 
analyzed archival X-Ray data for each of them. We have 
used Chandra and/or XMM data for all objects except 
A 520, for which only ROSAT data are available. 
We modeled the emission of all clusters in the 0.3 — 
7.0 keV band. Pixel values of all detected point sources 
were replaced with values interpolated from the sur- 
rounding background regions; the CI AO tools wavdetect 
and dmfilth were used for this purpose. 
Using the SHERPA software, we fitted the cluster surface 
brightness to elliptical 2-D /3-models (see Eq. [7} . Results 
are listed in Tabled Fitted models from Chandra and 
XMM observations are roughly consistent. 
The 25 clusters have a weighted median projected axis 
ratio of e pro j = 1.24±0.09, in very good agr eement with 
the va lue of (e pro j) = 1.25±0.19 obtained bv lMohr et al.1 
( 1995) from Einstein data of a lower-redshift sample of 
65 objects. Only six of the 25 clusters have a projection 
in the p.o.s. close to be circular (e pro j < 1.2). 

5.1. Circular Versus Elliptical (3 Models 

Although clusters are rarely circular in projection, 
some previous joint analyses of X-Ray and SZE data 
have assumed spherical symmetry. In order to bound the 
effects of this simplification, we have also modeled the 
surface brightness profiles of each sample cluster with a 
circular /3-model. The choice of circular rather than ellip- 
tical (3 model does not affect the resulting of the central 
surface brightness, as shown in the top panel in Fig. 
For a few clusters the fitted value of the slope /3 differs 
slightly between circular and elliptical models (middle 



panel of Fig. As would be expected, however, signif- 
icantly different values for the core radius are obtained 
with these two models (bottom panel of Fig. I B . This be- 
havior has already been noted bv lHuehes fc Birkinshawl 
(QUI. 

Therefore, relaxing the assumption of circular projection 
on the p.o.s. when measuring the angular diameter dis- 
tance fEa. 110(1 . mainly affects the value of the projected 
core radius # Cl proj- The bottom panel in Fig. ^ shows 
that the core radius obtained using a circular /3-model 
($c,Circ) is consistently lower (black squares) than that 
obtained from an elliptical model (# Ci eii)- $c,Circ can in 
fact be well approximated by the arithmetic mean of the 
two semi-axes of the elliptical isophotes in the p.o.s. 
The angular diameter distance obtained assuming spher- 
ical symmetry (Table |3J) can therefore, in first approxi- 
mation, be corrected for the observed ellipticity of the 
cluster in the p.o.s. multiplying -D c |exp by the correc- 
tion factor /coi- 

n iCirc.Corr _ n iCirc f _ n |Circ 1 + e proj /-, ,\ 
-^clExp — ^c| Exp /cor — -^clExp n • V 1 ^/ 

^proj 

As shown in the bottom panel of Fig. ^ the corrected 
values of the core radii c |cfrc 

^clcirc = ^c.Circ ~. (15) 

/cor 

provide a good approximation to C C \\ (gray squares). 

6. CLUSTER MORPHOLOGY IN THREE DIMENSIONS 

6.1. Angular Diameter Distances 

In order to estimate the l.o.s. extent of clusters, then, 
we need only to obtain values of -D c | Exp from the X-Ray 
and SZE data (via Ea. llO|) and compare them (via Ea. ll3fl 
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Table 2. Two-Dimensional Analysis 





















Cluster 


R. A. 


Decl. 






o 

(deg) 


r c 
(arcscc) 





Satellite 


MS 1137.5+6625 


11 40 22.3 


+66 08 15.3 


1 


113 ^ 0.014 


63.9 ± 1.0 


11.28 ± 0.55 


0.58 ± 0.01 


1 


MS 0451.6-0305 


04 54 11.4 


—03 00 51.3 


1 


.307 ± 0.015 


84.1 ± 1.1 


45.1 ± 1.2 


0.88 ± 0.02 


1 


CI 0016+1609 


00 18 33.5 


+ 16 26 12.9 


1 


.205 ± 0.013 


310.8 ± 1.7 


36.42 ± 0.93 


0.63 ± 0.01 


1 


CI 0016+1609 


00 18 33.1 


+ 16 26 10.5 


1 


.168 ± 0.019 


314 ± 3 


39.0 ±1.1 


0.65 ± 0.08 


2 


RXJ 1347. 5- 1145 


13 47 30.7 


— 11 45 09.1 


1 


.453 ± 0.019 


21.5 ± 1.0 


6.36 ± 0.15 


0.533 ± 0.003 


1 


A 370 


02 39 53.3 


— 01 34 39.0 


1 


.564 ± 0.018 


353.8 ±0.7 


50.1 ± 1.8 


0.52 ± 0.01 


1 


MS 1358.4+6245 


13 59 50.7 


+62 31 04.1 


1 


.325 ± 0.019 


23.4 ±1.4 


14.23 ± 0.49 


0.526 ± 0.004 


1 


A 1995 


14 52 57.9 


+58 02 55.8 


1 


.242 ± 0.010 


122.2 ± 1.0 


45.78 ± 0.88 


0.73 ± 0.01 


1 


A 611 


08 00 56.8 


+36 03 23.5 




1.14 ± 0.05 


326 ± 9 


21.84 ± 0.62 


0.596 ± 0.006 


1 


A 697 


08 42 57.6 


+36 21 56.8 


1 


.334 ± 0.016 


16.2 ± 1.2 


54.3 ±1.7 


0.64 ± 0.01 


1 


A 1835 


14 01 02.0 


+02 52 42.9 


1 


.225 ± 0.012 


7.0 ±1.4 


8.34 ± 0.14 


0.511 ± 0.002 


1 


A 2261 


17 22 27.1 


+32 07 57.4 


1 


.022 ± 0.017 


0.0 ± 1.7 


20.58 ± 0.75 


0.578 ± 0.007 


1 


A 773 


09 17 53.1 


+51 43 37.9 


1 


.237 ± 0.022 


0.0 ± 2.5 


49.5 ± 1.9 


0.63 ± 0.02 


1 


A 773 


09 17 52.7 


+51 43 37.0 


1 


.184 ± 0.019 


0.0 ± 2.9 


45.2 ± 2.3 


0.58 ± 0.06 


2 


A 2163 


16 15 46.6 


—06 08 44.9 


1 


.206 ± 0.004 


0.0 ± 0.6 


94.80 ± 0.95 


0.720 ± 0.005 


1 


A 2163 


16 15 46.0 


-06 08 46.9 


1 


.163 ± 0.009 


0.0 ± 1.4 


91.3 ± 1.3 


0.71 ± 0.01 


2 


A 520 


04 54 09.8 


+02 55 22.4 




1.06 ± 0.05 


347 ± 20 


115.2 ± 11.8 


0.59 ± 0.04 


3 


A 1689 


13 11 29.6 


-01 20 28.0 


1 


.141 ± 0.012 


342.4 ± 2.2 


27.01 ± 0.73 


0.582 ± 0.005 


1 


A 665 


08 30 57.1 


+65 51 01.8 


1 


.238 ± 0.012 


33.7 ± 1.2 


63.0 ± 1.3 


0.582 ± 0.006 


1 


A 2218 


16 35 51.9 


+66 12 34.6 


1 


.162 ± 0.009 


83.5 ± 1.5 


56.0 ± 0.84 


0.591 ± 0.004 


1 


A 2218 


16 35 52.4 


+66 12 33.5 


1 


.201 ± 0.013 


85.4 ± 1.7 


58.5 ± 1.1 


0.603 ± 0.009 


2 


A 1413 


11 55 17.9 


+23 24 16.2 


1 


.473 ± 0.019 


182.2 ± 0.9 


43.0 ± 1.4 


0.573 ± 0.007 


1 


A 2142 


15 58 20.1 


+27 14 03.5 


1 


.540 ± 0.007 


52.3 ± 0.3 


73.1 ± 1.1 


0.598 ± 0.005 


1 


A 478 


04 13 25.3 


+ 10 27 53.5 


1 


.477 ± 0.006 


43.5 ± 0.3 


30.53 ± 0.16 


0.533 ± 0.001 


1 


A 1651 


12 59 21.9 


-04 11 44.6 


1 


.184 ± 0.013 


272.4 ± 1.9 


73.4 ± 2.0 


0.598 ± 0.007 


1 


A 401 


02 58 57.1 


+ 13 34 37.8 


1 


.303 ± 0.008 


325.1 ± 0.6 


142.8 ± 3.0 


0.625 ± 0.007 


1 


A 399 


02 57 52.0 


+ 13 02 38.7 


1 


.207 + 0.009 


337.5 ± 1.2 


139.8 ± 3.8 


0.536 ± 0.008 


2 


A 2256 


17 04 00.4 


+78 38 37.1 


1 


.327 + 0.008 


61.3 ± 0.5 


529 ± 16 


1.33 ± 0.06 


2 


A 1656 


12 59 44.1 


+27 56 43.0 


1 


.141 ± 0.006 


0.0 ± 0.6 


540 ± 34 


0.65 ± 0.07 


2 



Note. — Fit parameters of the elliptical f3 model: x c ,y c is the central position; e pro j is the projected axial ratio; 9 is the orientation angle 
(north over east); r c and /3 are the model core radius and slope, respectively. In the last column, label 1 is for Chandra, 2 for XMM and 3 for 
ROSAT HRI observation. 



to the angular size distance obtained from the measured 
redshift and our adopted cosmological model. Since only 
the X-Ray data are publically available, however, we are 
unable to jointly fit both SZE and X-Ray data. For this 
reason we must rely on published values of central CMB 
temperature decrement (ATo) for our analysis. 
A potential difficulty with this approach i s that 
the available values of ATp, from iReese et alJ (2002) 
and iMason et alJ l|200if ). have been inferred assuming 
that clusters arc circularly symmetric when seen in pro- 
jection on the sky. While this assumption is quite rea- 
sonable given the limited spatial resolution of the data 
available to these authors, it is not, in general, consis- 
tent with the results of our analysis of higher-resolution 
X-Ray data (see TableEJ . In the limit of very high spatial 
resolution SZE data, we would expect this inconsistency 
to have negligible effect on our results, just as we find 
that with high-resolution X-Ray data, the same central 
X-Ray surface brightness is inferred from fits of circular 
and elliptical models (see the top panel in Fig. 
We have computed values of the angular diameter dis- 
tances for all clusters in the sample both under the as- 
sumption of spherical symmetry, and relaxing the as- 
sumption to a more general triaxial morphology. For 
the spherical case results obtained modeling the cluster 
X-Ray surface brightness profiles with circular /3-models 
(§ 15.1 \ were substituted into Eq. (|1C)|> . For the triax- 
ial case results from the elliptical /3-models were instead 
used (§EJ- The term f(v,T e ), which accounts for the 
frequency shift and also includes relativistic corrections, 



was computed as described bv lltoh et alJ l)1998fl . Central 
values of t he CMB temperatu re d ecrement (ATM wer e 
taken from IReese et all l|2002T) and IMason et al J (|200H) . 

The resulting values of -D c |ex P are listed in Table |3] The 
largest source of error (about 70% of the total) is the 
uncertainty on the SZE measurement of ATo. The sec- 
ond most significant error source is the uncertainty in the 
X-Ray measurement of the intra-cluster plasma temper- 
ature (about 20%). Both -D c |ex P ana - ^ c Iex P are plot- 
ted in Fig. [3 A comparison between the experimental 

quantities -Dc| E x P and -D c |ex P an d tlie values of -D c | C osm' 
together with their relative errors, highlights the high 
precision to which the cosmological angular distance is 
known, compared to the two experimental estimates, in 
support of our approach of a fixed cosmological model. 

6.2. Elongation Along the Line of Sight 
Assuming a general triaxial morphology, the ratio be- 
tween -D c | Exp an d ^clcosm provides an estimate of the 
ratio of the cluster axis along the l.o.s. and the cluster 
major axis in the p.o.s. We have computed values of 
-l7cIex P for all clusters in the sample f§ IB.lfl . For each 
cluster we have then also computed TJ> c | Cosm (Eq.0J) and 
have then estimated their ei. . s .- Resulting values are 
listed in Table H 

Since the observables in our analysis have asymmetric 
uncert ainties, we apply corrections given bv ID'Agostinil 
( 2004) to obtain estimates of sample mean and standard 
deviation. 
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Fig. 2. — ^-Mcosm as a function of redshift (solid line); dashed 
lines represent the ltr confidence level. Black and grey point rep- 
resent -D c |™ p ant ^ ^ c Iexp computed in this work, respectively. 
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Fig. 3. — Distribution of elongation along the l.o.s. for all clusters 
in our sample. 



Table 4. 3-D Morphology 



Table 3. Angular Diameter Distance 



Cluster 



Dc\ Co: 
(Mpc 



D, 



Exp 

(Mpc) 



d r" 

(Mpc) 



MS 1137.5+6625 
MS 0451.6-0305 
CI 0016+1609 
RXJ1347. 5-1145 
A 370 

MS 1358.4+6245 

A 1995 

A 611 

A 697 

A 1835 

A 2261 

A 773 

A 2163 

A 520 

A 1689 

A 665 

A 2218 

A 1413 

A 2142 

A 478 

A 1651 

A 401 

A 399 

A 2256 

A 1656 



1537±£ 
1322^76 
1318±7g 
1189lgg 
1063±g? 
9741" 
9641" 
893+f 
880+^ 

743+Ji 

799+30 
,zz -40 

686±H 

632 + ^ 
601±H 

3491^ 
S40±iJ 
327±» 
2&9t\ 2 6 
282±jg 
232±\° 
96+4 



3179 
1278 
2041 



,+ 1103 
-1640 
+ 265 
'-299 
+ 181 
'"'-■--514 
1991 +368 
lzzl -343 

M0Z -1245 

866 +if? 

1119_282 
995 + 325 

998+298 

iUZ '-198 

1049± 27 6 2 
1460±||i 

q n q + 18 1 

79 , + 270 
' zo — 236 

466+ 217 
1029±H» 

187l 2 7 2 

406 ±135 
3?3 + 202 

g 10 +693 

wtlt 

296tl 2 7 
235+^ s 



2479 ± 1023 
1073 ± 238 
1635 ± 391 
1166 ± 262 
1231 ± 441 
697 ± 183 
885 ± 207 
934 ± 331 
1099 ± 308 
946 ± 131 
1118 ± 283 
1465 ± 407 
806 ± 163 
387 ± 141 
604 ± 84 
451 ± 189 
809 ± 263 
478 ± 126 
335 ± 70 
448 ± 185 
749 ± 385 
369 ± 62 
165 ± 45 
242 ± 61 
103 ± 42 



Note. — Cosmologica l angular dia meter dista nce, its experi mental 
estimate as reported by IReese et alj 1200211 and [M a son et al.l 1200111 
assuming spherical symmetry, and the experimental quantity ^dj^ 1 , 
defined by Eq. J10I . computed in this paper assuming the clusters are 
oblate spheroids. 



All clusters in our sample were X-Ray selected; X-Ray 
surveys are surface brightness limited. Clusters close to 
the detection limit which are elongated along the l.o.s. 
will be detected, while the ones which are more extended 
in the p.o.s. will be missed. If a surface brightness limit 
is fixed which is far above the detection limit of the sur- 
vey, the prob lem should b e elimi nated. 
In bo th the IReese et ahl (|2JX)2f) and the iMason et all 
( 2001) samples this "correction" limit was applied. Our 
final sample shows in fact only mild signs of preferential 
elongation of the clusters along the l.o.s. (see Fig. |3J). 
Of the 25 clusters, 15 clusters are in fact more elon- 



Cluster 




ei.o.a. 






q„,ax 


MS 1137.5+6625 


1 


61 ± 


69 


1 


79 ± 0.68 


MS 0451.6-0305 





81 ± 


18 


1 


31 ± 0.11 


CI 0016+1609 


1 


24 ± 


30 


1 


49 ± 0.32 


RXJ1347. 5-1145 





98 ± 


22 


1 


45 ± 0.19 


A 370 


1 


16 ± 


42 


1 


81 ± 0.48 


MS 1358.4+6245 





72 ± 


19 


1 


40 ± 0.12 


A 1995 





92 ± 


21 


1 


24 ± 0.14 


A 611 


1 


05 ± 


37 


1 


19 ± 0.27 


A 697 


1 


25 ± 


35 


1 


67 ± 0.40 


A 1835 


1 


17 ± 


16 


1 


43 ± 0.18 


A 2261 


1 


51 ± 


38 


1 


54 ± 0.39 


A 773 


2 


03 ± 


56 


2 


51 ± 0.71 


A 2163 


1 


18 ± 


24 


1 


37 ± 0.23 


A 520 





56 ± 


20 


1 


77 ± 0.26 


A 1689 





95 ± 


13 


1 


14 ± 0.08 


A 665 





71 ± 


30 


1 


40 ± 0.29 


A 2218 


1 


85 ± 


44 


1 


57 ± 0.45 


A 1413 





93 ± 


25 


1 


47 ± 0.20 


A 2142 





96 ± 


20 


1 


54 ± 0.18 


A 478 


1 


32 ± 


54 


1 


95 ± 0.65 


A 1651 


2 


29 ± 1 


18 


2 


71 ± 1.39 


A 401 


1 


28 ± 


21 


1 


67 ± 0.27 


A 399 





58 ± 


16 


1 


71 ± 0.12 


A 2256 


1 


04 ± 


20 


1 


38 ± 0.22 


A 1656 


1 


08 ± 


43 


1 


23 ± 0.32 



Note. 
ratio. 



Cluster elongation along the l.o.s. and maximum axial 



gated along the l.o.s (ei. . s . > 1), while the remaining 
10 clusters are compressed. The mean of the distribu- 
tion of the elongations is (ei. . s .) = 1.15±0.08. In pres- 
ence of likely outliers, th e median is a more stable esti- 
mator ijGott et alJl200H) . The median of the ei. . s .'s is 
1.08 ± 0.17. While on average we observe only a very 
slight preferential elongation of the clusters along the 
l.o.s., residual of X-Ray selection effects, clusters with 
extreme axes ratios are still preferentially selected if the 
elongation lies along the l.o.s. This is a clear example of 
how deeply X-Ray selected cluster samples are affected 
by morphological and orientation issues. 

6.3. Maximum Axis Ratio 

We can estimate the three ellipsoidal axis lengths (v± , 
i>2 and i>3 ) from the measured values of ei. . s . arid e pro j, 
and from these the ratio of the semi-major to the semi- 
minor axis, g max - ? m ax is an extremely convenient tool 
to describe the intrinsic shape of a cluster since it allows, 
without the aid of further parameters, to quantify how 
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far a cluster is from spherical symmetry. For most clus- 
ters in the sample, the confidence regions of Vi, vi and 
i>3 are highly overlapping so that, for example, the upper 
bound of the \-a interval for the estimate of the ratio 
between the median and the minor axis, g m id, may be 
larger than the upper limit of g max - 
To obtain well defined estimates of the errors of the max- 
imum, intermediate and minimum axis ratios for each 
cluster, assuming the ?Vs to be normally distributed, we 
have obtained 10 5 random samples from each distribu- 
tion. We have then selected the maximum, the interme- 
diate and the minimum values of each set of three in or- 
der to build the distribution of the maximum, intermedi- 
ate and minimum axis ratios. We have finally computed 
the standard deviations of such three distributions, that 
provide estimates for the errors for the axes ratios. The 
resulting values of q m ax are listed in Table 01 and their 
distribution is shown in Fig. 0] 

9max has a mean value of (g lnax ) = 1.59±0.07, and a 
median of 1.49 ± 0.17. This result is consistent with 
cosmological simulations in which the me an value of 
the maximum axial ratio ra nges fro m 1.56 l|Suwa et al] 
120031 IKasun fc Evrardll200l to 1.8 ll.Ting fc Sutol l20021. 
The intermediate axis ratios, q m id show a median of 
1.21 ± 0.12. At 1 — a no cluster in the sample can be 
approximated as spherical; 11 clusters are spherical at 
the 3 — a confidence level. 

Although our estimates are affected by large errors and 
the data sample is of modest size, we look for trends in 
the distribution of the maximum axial ratios. 
No correlation is observed between the maximum axial 
ratio and redshift (see Fig. [S] where the solid and dashed 
lines represent the weighted and non weighted linear best 
fit to the data, respectively). 

A poor correlation is observed also between the maxi- 
mum axial ratio and the cluster gas temperature. We 
find at most a weak tendency for hotter clusters to ex- 
hibit smaller values of q ma x ■ The linear weighted best fit 
to the data is: g max = (2.02±0.35) - (0.060±0.041) T. 
The trend is plotted in Fig. El where the solid and dashed 
lines represent the weighted and non weighted linear best 
fit to the data, respectively. The absence of such a cor- 
relation may indicate that, in our sample, high cluster 
temperatures are not predominantly the r esult of shocks 
associ ated with accretion of sub-clusters Randal l et al.l 
(2002), since such accretion events seem likely also to 
produce departures from spherical morphology. 
From the distribution of axial shapes of clusters in our 
sample, we can estimate the effect that the assumption 
of spherical symmetry has on the determination of the 
total cluster mass. If the mass is computed at large dis- 
tances from the cluster center (> 1 Mpc) the difference 
between the two models is less than 2% even for the most 
elongated clusters. If the mass is computed close to the 
cluster core the effect becomes larger, ranging from 4% 
to 25%, for less to more elongated clusters in our sample, 
respectively, when the mass is computed within a sphere 
of radius 100 kpc. Triaxial cluster distributions could 
therefore at least partially account for the observed dis- 
crepancies in the total mass of clusters computed with 
lensing and X-Rays measurements. 

We then analyze a subsample of the 10 clusters for which 
the presence of a cooling flow has been claimed (i.e. for 
which the upper limit, 90% confidence, to the central 
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Fig. 4. — Distribution of maximum axial ratios. The gray his- 
togram is for cooling-flow systems. 
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Fig. 5. — Maximum axial ratio as a function of the redshift. The 
solid and dashed lines represent the linear best fit to the data with 
and without weights, respectively. 



cooling time has been measured to be less than 10 10 yr). 
Cooling flow clusters are typically recognized as dynam- 
ically relaxed systems in which the ICM is supported 
by thermal pressure which dominates over non-thermal 
processes. Their X-Ray emission is in most cases regular 
and symmetric and little or no substructures is visible 
at optical wavelengths. We find no indication, however, 
that cooling flow clusters are more likely to be spheri- 
cal. Fig. 01 suggests that the distribution of maximum 
axial ratios for the cooling flow sample is indistinguish- 
able from that of the sample as a whole; a Kolmogorov- 
Smirnov test confirms this impression. 
Finally, we find no relationship between cluster elonga- 
tion along the l.o.s. and 2-D ellipticity (see Fig. 01. In 
particular, a circular (projected) surface brightness pro- 
file is not an indicator that a cluster is in fact spherical. 

6.4. Ellipticity and Prolateness 

Triaxial ellipsoids can b e represented in the ellipticity- 
prola teness plane (E, P) ijThomas et al.lll998t iSulkaner] 
1999). The ellipticity is defined as: 



E 



e 2 - e 2 ■ 

L max ^min 



2e • + e 

mm 1 



2 

mcd 



and the prolateness as: 



i 2 „ 2(0 2 
p _ mm med 



^ ^min 



4- e 1 



— ^med — *m 



(16) 



(17) 



The 



where the axial ratios satisfy e n 
allowed region in the (E, P) plane is a triangle delim- 
ited by the lines on which prolate and oblate clusters fall 
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Fig. 6. — Maximum axial ratio as a function of temperature. 
The solid and dashed lines represent the linear best fit to the data 
with and without weights, respectively. 
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Fig. 7. — Cluster maximum axial ratio as function of (projected) 
cllipticity in the p.o.s. Grey squares denote cooling-flow clusters. 
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Fig. 8. — Distribution in the ellipticity-prolateness plane for our 
sample. Grey squares highlight cooling-flow clusters. 



(P = —E and P = E, respectively) and the line con- 
necting their endpoints. Fig. [5] shows the distribution in 
ellipticity-prolateness for our sample. No cluster in our 
sample shows extreme values of the elli pticity parameter. 
As e xpected from some simulations ijKasun fe Evrardl 
120041) . prolate shapes (P > 0) may be more likely (18 
clusters) than oblate ones (P < 0). Once again, cooling 
flow clusters (gray boxes) are indistinguishable from the 
sample as a whole. 

7. SUMMARY AND DISCUSSION 

In this paper we have discussed how observations of 
clusters in the microwave and X-Ray spectral bands can 
be combined to constrain their intrinsic 3-D shapes, pro- 
vided that the cosmological model is known. We have 
applied our method to a combined sample of 25 clusters 



of galaxies. In doing so we make the simplifying assump- 
tion that the clusters are ellipsoids with one axis parallel 
to the l.o. s. 

Our sample clusters were originally selected on the ba- 
sis of X-Ray luminosity, with the selection threshold well 
above the detection limit, in order to avoid selection on 
the basis of X-Ray surface brightness. Even so, we ob- 
serve that clusters with extreme axes ratios are still pref- 
erentially selected only if the elongation lies along the 
l.o. s. 

The mean value of the axial ratio in l ow-density cosmo- 
logical simulations ran ges from 1.56 HSuwa et al]l2003 
iKasym fc Evrardl l2004f I to 1.8 l|,Ting fc Sutol 120021) and 
2.0 ijThomas et al.lll998j) . This is consistent with results 
we present here (g max ) = 1.59±0.07. The spherical hy- 
pothesis is generally rejected, with prolate-like shapes be- 
ing slightly more likely than oblate-like ones. Numerical 
investigations suggest there should be some tendency for 
axial ratios to be large r at higher redshift llJing fc Sutol 
120021 iSuwa et al] 120031 iKasun fc Evrardl I2004J) . though 
this effect is only marginal: the axial ratio increase s only 
3% from z = to z = 0.8 l|Kasun fc Evrardll200l . Our 
data are not sufficiently precise to test this prediction. A 
poor correlation is also observed between the maximum 
axial ratio and the cluster gas temperature. The absence 
of such a correlation may indicate that high cluster tem- 
peratures are not mainly the result of shocks associated 
with accretion of sub-clusters, since such events would 
likely produce departures from spherical morphology. 
The uncertainties on our results are mainly due to the 
relatively large errors in the SZE measurements of cen- 
tral temperature decrement, together with the quadratic 
dependence of the distance estimates on this parameter. 
More accurate SZE measurements, with better effective 
angular resolution, are required to extend this analysis 
to a large sample spanning a greater redshift range. 
A relevant number of cosmological tests are today based 
on the knowledge of the mass of galaxy clusters through 
X-Ray measurements; these masses are though usually 
computed assuming a spherical symmetry. It is there- 
fore extremely important to assess the effect that such 
assumption has on the determination of the total cluster 
mass. We have estimated that while the effect is negli- 
gible when the mass is computed at large distances from 
the cluster center, the discrepancy between the two mass 
values becomes much more important as we get closer to 
the cluster core (i.e. ~ 25% when the mass is computed 
within a radius of 100 kpc). Triaxial cluster shapes may 
therefore at least partially account for the discrepancies 
between cluster mass computed with strong lensing and 
X-Rays data. 

Although the presence of a cooling flow is often inter- 
preted as a sign of dynamical relaxation, the properties 
of the subsample of 10 cooling flow clusters do not dif- 
fer from those of the whole sample: cooling flow clusters 
therefore do not show preferentially spherical morpholo- 
gies. 
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APPENDIX 

A. TRIAXIAL ELLIPSOIDS 

We consider a cluster electron density distribution described by an ellipsoidal triaxial /3- model. In a /3-model, 
the electron density of the intra-cluster gas is assumed to be constant on a family of similar, concentric, coaxial 
ellipsoids. High resolution TV-body simulations have shown the asphericity of density profiles of dark matter halos 
and how such profiles can be accurately described by concentric triaxial ellipsoids with aligned axis: ellipsoidal (3- 
models therefore provide a m ore detailed descrip tion of relatively relaxed simulated halos, respect to the conventional 
spherically symmetric model ijJing fc Sutoll2002() . 

In a coordinate system relative to the cluster, we then describe the cluster electron density as: 



-3/3/2 



n e = n e0 \1 + ±±i=i_p^ I (Al) 

where {a^mt}, with i = 1,2,3, define an intrinsic orthogonal coordinate system centered on the cluster's barycenter 
and whose coordinates are aligned with its principal axes; r c is the characteristic length scale of the distribution, which 
in our case is defined as the core radius; along each axis, Wj is the inverse of t he co rresponding core radius in units of 
r c ; n e Q is the central electron density. The electron density distribution in Eq. IjAlfl is described by 5 parameters: n e o, 
/3, the axial ratios ei = v\/vz and e2 = V2/V3, and the core radius r C 3 = r c /v^, along a^mt: 

(22 2 2 2 \ 3/3/2 
, . e l X l,int + e 2 x 2,int + ^.int \ / \ r>\ 

1 H -3-! '■ I (A2) 

To write the electron density distribution given by Eq. (|A2I) in a coordinate system relative to the observer, three 
additional parameters are needed: the rotation angles -#Eu> V'Eu and i/'Eu- of the three principal cluster axes respect 
to the observer. A rotation through the first two Euler angles is sufficient to align the a^obs-axis of the observer 
coordinates system {xi i0 b s } with the l.o.s. of the observer, i.e. the direction connecting the observer to the cluster 
center. When viewe d from an arbitrary direction, quantities constant on similar ellipsoids project themselves on similar 
ellipses l|Starkll977j) . A third rotation -i/'Eu- will align a^i j0 b s and a^.obs with the symmetry axes of the ellipses projected 
on the p.o.s. of the observer (p.o.s.). Eight independent parameters are therefore required to uniquely geometrically 
characterize the electron density distribution of a triaxial galaxy cluster. 

The axial ratio of the major to the minor axes of the observed projected isophotes, e pro j(> 1), is a function of the 
sha pe param eters and of the direction of the l.o.s., defined by the first two Euler angles, 6*Eu and tpEu- It is given 
by <|BinggeialT980|) : 



--proj 




J+1+ y/(j - 2 + 4fc 2 



V(j-0 2 +4fc 2 ' 



(A3) 



where j, k and I are: 



j = e\e\ sin 2 9 Eu + e\ cos 2 9eu cos 2 <^e u + e 2 cos 2 #e u sin 2 <^e u (A4) 

fc = (e 2 — e 2 ) sin^Eu cos (yCEu cos6*eu (A5) 

/ = e 2 sin 2 ip Eu + e\ cos 2 (^eu (A6) 

The rotation angle betwee n the principal axes of the observed ellipses and the projection onto the sky of the ellipsoid 
^3,int-axis is l)Binnevlll985f) : 



ib — — arctan 
2 



2fc 



(A7) 



U-i. 

The apparent principal axis that lies furthest from the projection onto the sky of the £3 ; i n t -ellipsoid axis is the apparent 
major axis if ijBinnevI 11985) 

(j - I) cos 2ip + 2k sin 2ip < (A8) 

or the apparent minor axis otherwise. In what follows, we assume a:i j0 b s to lie along the major axis of the isophotes, 
so that: 

2fc 



V'Eu = arctan 



j-i- ^fc 2 - (j - iy 



(A9) 



The projected axial ratio e pro j, the orientation angle of the projected ellipses, the slope (3 and the projection in the 
p.o.s. of the core radius can be determined fitting observed images to the /3-model; further independent constraints 
are therefore still required to uniquely determine the gas distribution. 

The X-Ray surface brightness and the SZ temperature decrement are given by projection along the l.o.s. of two 
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Fig. B9. — The error in the estimate of the axial ratio, g ma x, due to inclination issue vs. the triaxiality degree. The horizontal line is 
fixed to the true value of g max . Points with error bars refer to estimates in the hypothesis of one principal axis oriented along the l.o.s. 



different powers of the electron density n e . Foll owing Sta3 fll977j) . we calculate the projection along the l.o.s. of the 
electron density distribution, given by Eq. (£3, to a generic power m which, in the observer coordinate system, can 
be written as: 

^ ) ( 0) 

where D c is the angular diameter distance to the cluster and 9i = Xi jC ,bs/-Dc is the projected angular position on the 
p.o.s. of Xi^hs- h is a function of the cluster shape and orientation: 

h = e\ sin 2 #Eu sin 2 <^e u + e\ sin 2 #e u cos 2 i^eu + cos 2 $e u (AH) 

The observed cluster angular core radius # c ,proj is the projection on the p.o.s. of the cluster angular intrinsic core 
radius 

"e -^ 1/2 

Cproj 



?c3 



eie 2 



h 1 ^ (A12) 



where 6> c3 = r c3 /D c . 

B. INCLINATION ISSUES 

The analysis presented above is based on the main assumption that one cluster principal axis is elongated along 
the line of sight. Despite this assumption is quite strong, it does not affect significantly the results. To test the 
effect of inclination on the estimate of the axis ratios, we proceed in the following way. First, we generate a galaxy 
cluster, characterized by the maximum axis ratio q ma x and by a parameter related to the degree of triaxiality, T = 
(e m id — e nml ) / (e max — e m ; n ); oblate and prolate clusters correspond to T = and 1, respectively. Since we are only 
interested in inclination issues, we neglect other measurements errors. Then, we generate a set of 25 viewing angles 
{(?Euj VEu}- We assume that orientations are completely random, i.e. the angle #Eu is between and tt and follows 
the distribution sin#Eu/2, whereas tpEu follows a uniform distribution between and 2tt. For each pair of orientation 
angles, we compute the projected ellipticity e pro j and the elongation ei. . s . and, then, we obtain an estimate of the 
axial ratios of the cluster in the hypothesis of one principal axis being aligned along the l.o.s. Finally, we calculate 
mean and standard deviation of the set of axial ratios corresponding to different viewing angles. If the value of such 
a mean is near that of the simulated cluster, then inclination issues hardly affect our analysis. In Fig. IJJ9I we plot the 
effect on the estimate of g max for different values of T in the case of q m & x = 1.5. Error bars equal standard deviations. 
As we can see, the error is minimum for highly triaxial clusters (T ~ 0.5), being ~ 0.1. Such an error should be added 
in quadrature to the value of A(7 max estimated in the previous sections, but due to its smallness it does not contribute 

significantly. The error increases for oblate or prolate clusters, with values ~ 0.4. In this paper we have been facing 
with the hypothesis of triaxial clusters, in fact the algorithm illustrated in the previous sections is optimized to describe 
triaxial spheroids. The capability of ellipsoids of revolution to reproduce the observed data set will be the subject of 
a forthcoming paper. This considerations are quite general and still holds for very different values of g ma x- The error 
due to inclination issue on the estimate of axial ratios can be usually neglected with respect to other uncertainties, in 
particular in the measurement of the SZ temperature decrement. 

C. GRAVITATIONAL LENSING 

Clusters of galaxies act as lenses deflecting light rays from background galaxies. In contrast to SZE and X-Ray 
emission, gravitational lensing does not probe directly the ICM distribution but maps the cluster total mass. The 
ICM distribution in clusters of galaxies traces the gravitational potential. Since we are considering a triaxial /?- 
model for the gas distribution, the gravitational potential turns out to be constant on a family of similar, concentric, 



12 



De Filippis et al. 



coaxial ellipsoids. Ellip soidal potentials are widely used in gravitational lensing analyses to fit multiple image systems 
(jSchneider et alJH992|) . 

The distribution of the cluster total mass can be inferred from its gas distribution. If the intra-cluster gas is assumed 
to be isothermal and in hydrostatic equilibrium in the cluster gravitational potential, while non-thermal processes are 
assumed not to contribute significantly to the gas pressure, the total dynamical mass density reads: 

*— - (sasb) y2 """•» (C1) 

where G is the gravitational constant and [im v is the mean particle mass of the gas. If we assume that the electron 
density of the ICM follows a /3-model distribution given by Eq. I)A2|I . the total gravitating mass density, in the 
coordinate system relative to the cluster, is given by: 



3f3k B T c 

Pai 



47rG^r7i p r^3 



'c3 



1=1 



' c3 1 + ^ellAcS 



(C2) 



where = 1 and r e ii is the ellipsoidal radius, r 2 u = ^j=i i e i x i,int) 2 - The projected surface mass density can 
subsequently be written, in the observer reference frame, as: 

-3/2 

^ = ( ' , 1 \ I 1 - ~TT~^ ! (C3) 




where: 

^_3 (3k B T cy ^l + e 2 proi 1 1_ 

4 Gprrip h 3 / 4 ^/e proj c ,proj D c 

Although the hypotheses of hydrostatic equilibrium and isothermal gas are very strong, total mass densities obtained 
under such assumptions can yield accurate estimates even in dynamically active clusters with irregular X-Ray mor- 
phologies. Elliptical potentials motivated by X-Ray observations were employed in the irregular cluster AC 114 to 
provide good fit to multiple image systems ifDe Filippis et alJl2004fl . 
The lensing effect is determined by the convergence k: 

k = (C5) 



which is the cluster surface mass density in units of the surface critical density E c 

c 2 A 



AttG D c D c 



(C6) 



where D cs is the angular diameter distance from the lens to the source and D s and D c are the angular diameter 
distances from the observer to the source and to the lens, respectively. Fitting the observed surface mass density to a 
multiple image system, it is possible to determine the central value of the convergence: 



which, using Eqs. (|C4|) and (|C6fl . can be written as 



koiZs) = 3^/3fc B T c yHT^l + e 2 roj l D cs 
<?[Lm v h 3 / 4 ^/eproj 6> c , proj D s 
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